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Field emission (FE) and the electronic-states origin of graphitic nanostructures were investigated 
by first-principles calculations based on time-dependent density-functional theory. We find that 
the FE current from graphitic ribbons changes remarkably depending on the hydrogen termination 
and the direction of the applied electric field. Also, the FE current from graphene sheets shows 
a dramatic increase around vacancy defects. We verified, through the analysis of local electronic 
structures and energy distributions of emitted electrons, that the dangling-bond (or a) character is 
responsible for these results and governs the nature of the FE of graphitic nanostructures. 

PACS numbers: 79.70.+q, 61.48.+C, 71.15.Mb 



I. INTRODUCTION 

Recent advances in nanotechnology have made possible 
the fabrication of a wide variety of nanometer-scale de- 
vices using graphitic nanostructures, such as fullereneSfi 
carbon nanotubes^ and other carbon-based materials^*^ 
Among the various applications for these nanostructures, 
electron field emitters made from carbon allotropes show 
significant promise for novel electronic devices f£*£***£ be- 
cause they can maintain stable forms under extremely 
high field-emission (FE) current densities owing to tight 
covalent-bonds. Chen et al. performed an FE experiment 
using graphite platelet nanofibers (GPNs)* 9 - with several 
thousand graphitic ribbons stacked together like a deck of 
cards. 3 Wu et al. carried out FE measurements of carbon 
nanowalls (CNWs), which are nano-graphitized sheets 
grown perpendicularly on substrates^** 0, and pointed out 
that CNWs are good candidates for nanoscale field emit- 
ters as well as carbon nanotubes— For the development 
of efficient field emitters, examination of FE properties 
of graphitic nanostructures is interesting and meaningful 
both from the viewpoint of fundamental science as well 
as technological applications of the field emitters. 

Over the past few years, interesting phenomena in 
graphitic ribbons have been predicted by theoretical 
studies— i* 2 **^** 1 - 4 ***^ According to theoretical results based 
on the tight-binding model for the tt electrons jiLi2iii, 
zigzag graphitic ribbons are always metallic, while arm- 
chair graphitic ribbons arc either metallic or insulating 
depending on the widths of the ribbons. The study also 
pointed out that zigzag ribbons exhibit specific electronic 
states, highly localized at the ribbon edges (edge states). 
This was supported by first-principles calculations 1 - 4 ***^ 
based on density- functional theory (DFT) — For the edge 
state, magneticfii* 1 * 3 " optical^ and thermal^ properties 
of graphitic ribbons have also been studied theoretically. 
However, thus far, there have been very few experimental 
studies of graphitic ribbons. 19-20 

Recently, the FE of graphitic ribbons was investi- 
gated by Tada and Watanabe using the time-dependent 



density-functional theory (TD-DFT) .21 They found that 
the dangling-bond states and not the edge states con- 
tribute primarily to the FE. They also found that not 
only evaluation of work functions but also knowledge of 
local electronic properties, the a- or 7r-bonding states, 
are prerequisite for understanding the microscopic mech- 
anisms of FE properties of covalently-bonded nanos- 
tructures. Tada and Watanabe's results cannot be de- 
rived from the conventional Fowler-Nordheim theory 2 *^ 
for the FE from flat jellium surfaces with the free- 
electron approximation. With the ultimate goal of com- 
plete understanding of microscopic mechanisms of the 
FE of graphitic nanostructures, in this study, we have 
investigated the effects of the electric-field direction and 
vacancy defects on the FE characteristics of graphitic 
nanostructures by using TD-DFT calculations. 

The outline of this paper is as follows. In Sec. II, the 
computational methods and models used are described. 
In Sec. Ill, we present numerical results for the effects 
of hydrogen (H) termination and the direction of applied 
electric field on the FE of graphitic ribbons. We also 
present results for the effect of vacancy defects on the 
FE of graphene sheets. We conclude with Sec. IV. 



II. METHODS AND MODELS 

We carried out conventional DFTi& and TD-DFT- 2 ^ 
calculations to investigate the electronic states and 
electron emission properties of graphitic ribbons and 
graphene sheets. TD-DFT— has been successfully 
applied to various non-equilibrium electron dynamics 
phenomena* 2 - 4 * 2 *^ including FE of carbon nanotubes* 2 *! 
graphitic ribbons* 21 - covalently-bonded atomic wires* 2 *^* 
and diamond surfaces p 2 *!*?*! We employed the norm- 
conversing pseudopotentials of NCPS97 31 « based on the 
Troullicr-Martins algorithm 32 and the generalized gradi- 
ent approximation by Perdew et al. for the exchange- 
correlation potential— The block Davidson algorithm 
was adopted to diagonalize the Kohn-Sham Hamiltonian 
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matrix. The electronic wave functions were expanded in 
terms of a plane-wave basis set. In the TD-DFT calcula- 
tions, the expansion coefficients become time-dependent 
quantities. 

First, we determined the ground states of atomic and 
electronic structures by using conventional DFT in a zero 
electric field. Next, we applied an electric field of 10 
V/nm and calculated the time evolution of the wave func- 
tions by applying the TD-DFT scheme. Then, wave func- 
tions with higher energies start to contribute to FE. The 
KS Hamiltonian is updated by a new electron density 
at each time step, and a part of the dynamical screen- 
ing effect and some electron correlation effects are au- 
tomatically taken into account in this calculations. It 
is noted that the simulation time should not exceed a 
certain time because the emitted electrons bounce back 
from the unit-cell boundaries. However, within a simu- 
lation time t < 150 a.u. (3.6 fs), such reflection does not 
occur. In this calculation of the time evolution of wave 
functions, we have adopted the seventh-order Taylor ex- 
pansion method. 2 A The Taylor expansion method cannot 
guarantee the unitarity condition, although this method 
has computational simplicity. The Suzuki- Trotter type 
split-operator method^* would be more accurate and re- 
liable for long time simulations^ because this method 
ensures the unitarity condition for the time evolution. 
The Taylor expansion method, however, guarantees the 
normalization condition for the electron number with an 
accuracy of 10 -6 for our simulation time. Moreover, the 
computational time for the Taylor expansion method is 
much shorter than that for the Suzuki- Trotter expansion 
method. Therefore, we have employed the Taylor ex- 
pansion method in this study. Finally, we counted the 
number of electrons tunneling into a large vacuum re- 
gion by summing the squares of the coefficients at each 
time step and obtained a curve of the electron number 
as a function of time. The value of the FE current was 
evaluated from the linear slope in the curve. 

The systems we investigated in this study are zigzag 
graphitic ribbons with and without H termination and 
graphene sheets with and without a vacancy defect, as 
shown in Fig.^ Since the effect of edge states on the FE 
is interesting, we chose the zigzag ribbons in the present 
calculations. The x and y axes are defined as shown in 
Fig. ^ The unit cell sizes for the graphitic ribbons are 
79.8 A x 2.49 A x 4.99 A and 14.8 A x 2.49 A x 88.6 
A, when electric fields are applied parallel (x) and per- 
pendicular (z) to the ribbons, respectively. The widths 
of the graphitic ribbon in Fig. Q^a) and (b) are 9.12 A 
and 7.10 A, respectively. The size of the unit cell for 
the graphene sheets is 8.52 A x 7.38 A x 85.2 A. The 
cutoff energies for the calculations for graphitic ribbons 
and graphene sheets are 45 Ry and 34 Ry, respectively. 
We chose 10 k y points in the first Brillouin zone for the 
graphitic ribbons and 16 points in the k x -k y plane in 
the first Brillouin zone for the graphene sheets for the 
determination of the ground-state electronic and atomic 
structures. The atomic structure around the vacancy de- 



fect of the graphene sheet was noticeably deformed, as 
seen in Fig.^d). 



III. RESULTS AND DISCUSSIONS 
A. Direction of electric field 

We begin with a brief review of a prior studySi on the 
FE of zigzag ribbons with and without H termination 
in a parallel electric field, £y/, along the x direction in 
Fig. ^ This previous study clarified the effects of H- 
termination on the FE properties of the zigzag ribbons 
using local electronic structures. The results of energy 
distributions of the emission currents and the electronic 
states responsible for the peaks are shown by broken lines 
in Fig. |3 For the H-terminated ribbon (Fig. 13a)), there 
are two peaks originating from a states in the energy 
distribution. The electronic distribution of emitted elec- 
trons resulting in the higher peak is shown in the up- 
per right panel in Fig. Ufa). The edge state (tt state) 
is not a contributor to the FE of the H-terminated rib- 
bon, even though the edge state is near the Fermi level 
for the H-terminated ribbon. On the other hand, for 
the clean ribbon (Fig. [2Jb)), a prominent peak appears 
in the energy distribution. The electronic-states origin 
of the sharp peak can be gleaned from the upper right 
panel in Fig. Gib) to be the dangling-bond (DB) state. 
The reason why electrons are emitted from the a state 
(H-terminated) or the DB state (clean) , and not from the 
7r states, can be elucidated from the directional angle be- 
tween the local electronic distributions and the applied 
field, Eij . Typical electronic orbitals, tt, a and DB around 
an edge of the clean ribbon are schematically shown in 
Fig. 13 As seen in Fig. |3 a or DB are parallel to the 
parallel field, E//, and tend to respond to the field. 

To verify that the electronic states contribute to the FE 
as long as the electronic orbitals protrude along the di- 
rection of the electric field, we have investigated whether 
7r orbitals react to an electric field (E±, z axis in Fig.^) 
that is perpendicular to the ribbon sheet. 

We now describe the results for the FE of the zigzag 
ribbon in the perpendicular field, E±. The energy distri- 
butions of the FE current and the corresponding elec- 
tronic distributions are given by solid lines in Fig. |21 
For the H-terminated ribbon (Fig.^a)), FE currents are 
emitted from the tt states at the edge, as expected, al- 
though no prominent peaks appear, compared to the en- 
ergy distribution curve for the parallel field, Eii (broken 
line in Fig.^a)). For the clean ribbon (Fig. Elb)), how- 
ever, a sharp peak originating from the DB state appears, 
similar to the case for the parallel field, Eu. Negligibly 
small contributions of the 7r states to the total FE cur- 
rent are found in the energy distributions. The direction 
of electronic distribution of the DB state in the upper 
right panel in Fig. |2b) is not parallel to the perpendic- 
ular field, E±. It follows from this result that the DB 
state contributes to the FE even when the direction of 



3 



the electronic distribution is not parallel to the direction 
of the electric field. 

Here, it is important to clarify why 7r states cannot be 
the main source of the FE even when it orbitals extend 
along the applied electric field, E± (See Fig. 0). Fig. 0] 
shows the equipotential surface of the clean ribbon on 
the plane perpendicular to the ribbon sheet (inset) and 
electrostatic potential energy curves along the direction 
of the perpendicular field, E± . It becomes apparent from 
the equipotential surface that the electric field becomes 
negligibly small due to the screening effect in the center 
of the ribbon. Therefore, electrons are not emitted from 
the center region of the ribbon. On the other hand, the 
electric field is enhanced at the edges of the ribbon, as 
seen from the equipotential surface, and thus electrons 
are emitted mainly from the edges. We focus on the 
three electrostatic potential energy curves drawn along 
the lines corresponding to the three positions in the inset 
of Fig. 0] The slope of the potential curve around the 
origin of the horizontal axis is the strength of the local 
electric field around the edge atoms along the direction 
of the perpendicular field, E± . The direction of the local 
field of curve 1, where the DB orbital exists, is parallel to 
the applied field, E±. However, the direction of the field 
for curve 3, where ir orbitals protrude, becomes oppo- 
site to that of the applied field. Consequently, electrons 
cannot be emitted from the n orbitals even when the it 
orbitals are parallel to the applied field. The electrostatic 
potential energy curves of the H-terminated ribbon also 
show similar features as the clean ribbon. 

Finally, we can summarize the FE properties of zigzag 
graphitic ribbons. We list the values of the FE currents 
from the two type of ribbons, H-termination and clean, 
in electric fields along two directions [E/i or Ej_) in Ta- 
ble HJ The DB orbital is the main source of the FE from 
the clean ribbon in both electric field directions. The DB 
state disappears upon H-termination and thus the a or- 
bitals and the 7r orbitals become emitting sources in the 
parallel and perpendicular electric fields, although the 
contributions of the 7r orbitals are negligible. 

B. Vacancy defects 

Having clarified the important role of the DB state in 
the FE from the graphitic ribbons, we further explored 
the FE property of the graphene sheet with an atomic 
vacancy (Fig.QJd)), because the DB states are generated 
in the vacancy. The graphene sheets investigated are 



TABLE I: FE currents pA/(unit cell) from two type of rib- 
bons (H-termination or clean) in electric fields (10 V/nm) 
along two directions {En or E±). 



Type of termination 




E± 


Clean 


0.37 


0.10 


H-termination 


0.14 


0.06 



shown in Fig.^c) and (d). The electric field was applied 
in a direction perpendicular (E±) to the sheet. 

We investigated the FE current-density images of the 
graphene sheets to clarify the role of the vacancy on FE 
characteristics. Figure |S] shows the FE images for the 
graphene sheet (a) without and (b) with vacancy defects. 
The plane on which the FE current distribution is plot- 
ted is 4.9 A above the graphene sheet. The FE current 
density around the vacancy is about 22 times as large as 
that from the other regions. The total FE currents ob- 
tained from the graphene sheet without and with vacancy 
defects are 0.054 and 0.106 //A/(unit cell), respectively. 
The FE images thus obtained reflect a dramatic change 
in the electronic distribution around the vacancy defect 
due to the removal of a carbon atom. It becomes appar- 
ent from the FE images that vacancy defects significantly 
enhance the FE current density of the graphene sheet due 
to the appearance of the DB states. It should be noted 
that the FE image in Fig. |3Jb), which is taken 4.9 A 
above the graphene sheet, does not necessarily indicate 
atomic-scale resolution FEM because the electron beams 
from each vacancy site interfere with each other before 
reaching the screen that is far from the emitter. 



IV. CONCLUSIONS 

The FE and the electronic-states origins of graphitic 
nanostructures were investigated by first-principles cal- 
culations based on TD-DFT. We found that the character 
of the local electronic states responsible for FE changes 
and thus the FE current varies, depending on the con- 
ditions of hydrogen termination and the direction of the 
electric field. The DB states are predominant sources for 
the FE because the electronic orbitals tend to appear at 
the edges and protrude along the direction of the elec- 
tric field in vacuum. High FE currents from the edges 
of the graphene sheets of CNWsjiS which are in paral- 
lel to the electric field, have been observed. Thus, the 
present study can provide a theoretical interpretation for 
the features observed in the experiment. We naturally 
expect that a similar property will be observed in the 
FE of GPNsiS It is interesting to note that the direction 
of electronic-orbitals has been found to remarkably influ- 
ence the ionization rate of atoms in strong laser fields^ 

Substantially enhanced emission current from the 
atomic vacancy sites in the graphene sheet confirms the 
essential contribution of the DB states to the FE from 
graphitic nanostructures. The 7r states, however, con- 
tribute minimally to the FE even when the electronic 
orbitals project along the direction of the electric field 
because the electric-field strength tends to decrease sub- 
stantially in the region of the 7r orbitals, i.e., owing to a 
screening effect. 

These findings are remarkable reflecting the covalent- 
bond character. The results emphasize the need for first- 
principles studies taking into account the geometric and 
electronic structures of field emitters towards better un- 
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derstanding of FE mechanisms, and for the design of 
nanoscale graphitic field emitters. 
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FIG. 1: Top views of graphitic ribbons (a) with and (b) 
without hydrogen (H) termination and graphene sheets (c) 
without and (d) with a vacancy defect in the unit cell. Gray 
circles are carbon atoms and the small black dots in (a) are 
H atoms. Dashed straight lines denote boundaries of the unit 
cells and a dashed circle in (d) denotes a vacancy defect. 

FIG. 2: Energy distribution of the FE current for (a) the 
H-terminated ribbon and (b) the clean ribbon. The vacuum 
level is chosen for the origin of the energy. On the right side 
of each panel, electronic distributions causing the peaks are 
shown by blue clouds with carbon atoms (white spheres) and 
H atoms (white dots) . The broken and solid lines denote the 
results in a parallel (E//) and perpendicular (E±) electric field, 
respectively. 

FIG. 3: Schematic diagram of electron orbitals (gray clouds) 
of a clean graphitic ribbon around the right edge (side view). 
Eii and E± denote the applied electric fields parallel and 
perpendicular to the ribbon, respectively. DB represents a 
dangling-bond. 

FIG. 4: The electrostatic potential energy as a function of 
the position from the ribbon sheet. The plane on which the 
ribbon lies is chosen to be the origin of the horizontal axis. 
The inset shows the equipotential surface of the clean ribbon 
from a lateral view in the perpendicular field, E±. Electrons 
are emitted from the ribbon into the upper region. The lines 
1, 2 and 3 are 1.2 A, 0.6 A and A away from the right edge 
atom of the ribbon. 

FIG. 5: FE images of the graphene (a) without and (b) with 
a vacancy defect. The distance between the FE image plane 
and the graphene is 4.9 A. The FE current density around 
the vacancy is about 22 times larger than that from the other 
regions. 
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